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Abstract. We reconsider non-thermal production of WIMP dark matter in a systematic way 
and using a numerical code for accurate computations of dark matter relic densities. Candidates 
with large pair annihilation rates are favored, suggesting a connection with the anomalies in the 
lepton cosmic-ray ux detected by Pamela and Fermi. Focussing on supersymmetric models we 
will consider the impact of non-thermal production on the preferred mass scale for dark matter 
neutralinos. We have also developed a new formalism to solve the Boltzmanns equation for a 
system of coannihilating species without assuming kinetic equilibrium and applied it to the case 
of pure Winos. 



1. General framework 

In this work we present a numerical code, now added to the package DARKSUSY [Ij, aimed 
to determine the dark matter rehc density when this is produced by the decay of of heavy 
long lived states with weak enough interactions with ordinary matter to be out-of-equilibrium 
also in the early stages of the history of the Universe. These states, referred for simplicity as 
cosmological moduli Xi, dominate the early stages of the history of Universe and lately decay 
into DM as well as SM particles with a sharp increase of entropy density. When this decay 
occurs at lower temperatures than the standard freeze-out temperature Tt,f,o. this non-thermal 
component dominates the dark matter relic density being the component already present in 
the thermal bath diluted by the entropy injection. Starting from the simplest case of only one 
decaying field without referring a definite particle physics scenario we have numerically recovered 
an interesting approximate trend occuring when, during the production process the dark matter 
pair annihilation rate F = ny^{av) is larger than the expansion rate H. In such a case pair 
annihilations are very efficient and instantaneously decrease the number density of x to the 
critical density level corresponding to F ~ when the annihilations stop. The non-thermal 
relic density is given by a simple scaling law as a function of the one obtained by assuming the 
DM candidate a thermal relic freezing-out at Tt,f_o.- 

with Trh is a characteristic temperature, numerically determined by the code, identifying the 
end of the decay process of the cosmological moduli and that can be analytically estimated 
by the condition F = H{Tiih) where F is the decay rate of the heavy field. Notice that a 
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Figure 1. Upper panel: Models fitting the cosmological value of the DM relic density in 
the plane {m^, {crv)). Colored regions are Pamela and Fermi best fit. The lines WMAP and 
Planck represent the current bound on the annihilation cross section from CMB and the Planck 
projected sensitivity. Lower panels: Curves of the cosmological relic density in two-dimensional 
slices of the Mi, M2, fi parameter space, for a few values of Trh- The filled areas correspond to 
regions violating the LEP lower bound on the chargino mass. 

DM candidate x whose thermal relic density is small compared to the DM density because the 
annihilation rate is too large, may become a viable dark matter candidate for an appropriate 
value of Trh- This suggest an intresting connection with anomalies detected by Pamela and 
Fermi which could be explained by an annihilating cross section 2-3 orders of magnitude larger 
than the thermal one. 

2. Application to the MSSM 

Let's now consider a definite particle physics scenario, i.e. the MSSM with neutralino Dark 
Matter. In the aim of determining the mass DM mass scale favored by non-thermal production 
we restrict to a class of models referred as split supersymmetry [2], characterized by very high 
mass scale for the scalar superpartners. In this way we could compute the non-thermal relic 
density for several values of the reheating temperature as a function of the Bino and the Wino 
mass parameters Mi and M2, and of the Higgs superfield parameter fi. 

The scenario depicted in the previous section is fitted by Wino-like or Higgsino-like dark 
matter. As reported of in fig. ([1]) the cosmological value of the DM relic density can be achieved 
for masses of the order of one hundred GeV, contrary to the thermal case where Wino and 
Higgsino masses in the TeV range are required. 



3. Kinetic equilibrium and decoupling 

The results presented until now are obtained under the assumption that DM is in kinetic 
equilibrium with the SM states of the thermal bath along all its production process. In case 
that DM is part of a system of coannihilating particles, this allows to describe it by a unique 
Boltzmann equation for sum of their number densities with an effective annihilation cross section 
defined as a weighted sum of the thermally averaged cross sections of the coannihilating particles. 
The validity of the kinetic equilibrium assumption is typically guaranteed by invoking a crossing 
symmetry among annihilations and scattering with SM particles. This assumption may not 
hold in all the models. Furthermore the picture is more complicated in presence of non- 
thermal production. First of all DM produced out-of-equilibrium at relativistic energies and 
thermalization with the heat but may not occur; secondly non-thermal freeze-out may occur at 
low temperatures, where the scatterings are less efficient. 

How to address this issue varies from one model to another. We have considered the case 
of a class of models referred as G2-MSSM p]. Here DM is represented by a system of two 
coannihilating species: the lightest neutralino and chargino (both wino-like) which are quasi 
degenerate in mass. The remaining spectum is at a very high mass scale, 0(10 — 100)TeV with 
the exception of the other gauginos, below the TeV scale, but non relevant for coannihilations. 
The issue of relaxing the condition of kinetic equilibrium have been approached in two steps. 
First of all we have verified the thermalization of the two species soon after production. This 
depends on their energy loss rates due interactions with the SM particles of the thermal bath. 
Thermalization occurs essentially instantaneously due to the following facts: at temperatures 
higher than order of 1 GeV charginos are promptly converted into neutralinos with the latters 
losing energy very efficiently. At lower temperatures charginos are able to lose energy before 
being converted and guarantee the thermalization of the two species system as long as inelastic 
neutralino scatterings are efficient. This picture is summarized in fig. ([2]). Once thermalization is 





Figure 2. Ratios between the effective rate of energy loss rate —1/E- dE/dt (dashed lines), 
or of the scattering/decay rate F (solid lines), to the Universe expansion rate for a few 
processes involving charginos and neutralinos. The panels refer to a G2-MSSM DM model with 
= 103.5 GeV and Trh = 100 MeV; the plot on the left hand-side refers to non-relativistic 
particles, E/m^ = 1.005, while that on the right-hand side corresponds to a sample relativistic 
case, E/rriy^ = 10. 

accomplished, the evolution of the system until kinetic decoupling have been followed by writing 
two Boltzmann equations for the distribution functions of the two species. These equation can 
be integrated by expanding the collisonal operator in an analogous manner to [1]: 
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with r^o^^i and r^±, being, respectively the scattering and the decay rate. The last term in the 
second equation is the non-thermal production term and is defined in [5j. The same formalism 
describes the dark matter distribution close to kinetic decoupling. Indeed, integrating the second 
moment of the distribution function it is possible to define the effective temperature T^q (we 
define only one temperature being charginos always kept in equilibrium by electromagnetic 
interactions) which coincides with the thermal bath temperature until there is kinetic equilibrium 
while scales as T^q oc after kinetic decoupling. In the second panel of fig ([3]) we show the 
behaviour of Ty.^ for different values of the reheating temperature. The impact of the decay of 
the cosmological moduli is twofold (fig. ([3])): the chargino decays tend to populate the system 
with neutralinos that are on average more energetic than for a thermal distribution, delaying 
the onset of the regime T^o oc and making the transition into this regime to be less sharp; at 
the same time, if Trh is so low that reheating increases significantly the expansion rate of the 
Universe H at the time of kinetic decoupling the departure from T^^ = T tends to be anticipated. 
In conclusion we point out that the developed formalism for the study of dark matter close to 
kinetic decoupling can be applied to other particle physics frameworks, also within standard 
cosmology. 




Figure 3. Left panel: ratio of the chargino number density over the neutralino number density 
for several values of Thh ■ Right panel: Ratio /T as function of the temperature of the thermal 
bath for same values of TptH ■ Plots are obtained for a Wino with mass equal to 200 GeV, however 
results depend only on the chargino-neutralino mass splitting which is about 160 MeV in the 
scenario under consideration. 
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